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Abstract 

We investigate transverse electromagnetic waves propagating in a plasma in the de Sitter space. 
Using the 3 + 1 formalism we derive the relativistic two-fluid equations to take account of the effects 
due to the horizon and describe the set of simultaneous linear equations for the perturbations. We use a 
local approximation to investigate the one-dimensional radial propagation of Alfven and high frequency 
electromagnetic waves and solve the dispersion relation for these waves numerically. 
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1 Introduction 



In recent years there have been renewed interests in investigating plasmas in curved spacetimes of general 
relativity; because, a successful study of the waves and emissions from plasmas falling into a compact body 
(e.g. black hole) will be of great value in aiding the observational identification of black hole candidates. 
The de Sitter (dS) space with positive cosmological constant has properties similar to a black hole. We 
study two-fluid plasmas near the horizon of the pure dS space. 

Over the last few decades, physicists have a growing interest in dS space. In the 1970s, the attention 
was due to the large symmetry group of dS space, which made the field theory in dS space less ambiguous 
than, for example, in the Schwarzschild spacetime. Researches in the 1980s focused the role it played 
during inflation-accelerated expansion in the very early universe. The universe is currently asymptotic 
dS and approach a pure dS space. Recent cosmological observations PQ, [21 [31 SI 0, E] suggest the possibility 
of existing a positive cosmological constant (A > 0) in our universe and this possibility gives the picture, 
among many others, of some features closely related to black holes: the existence of cosmological event 
horizons. These causal horizons exist even in the absence of matter, namely in empty dS space, and hide 
all the events which are not accessible for geodesic observers. In addition, the success of the ADS/CFT 
correspondence [8] [9] [10] has led to the intense study of dS space in the context of the quantum 
gravity [11]. The attention has been to obtain an analogue of the ADS/CFT correspondence in dS 
space, i.e. dS/CFT correspondence Q]3[Il[Il[IIE[]Il[rn[]13[I^ in the 

light of which there has been an extensive study of the semiclassical aspects of dS and asymptotic dS 
spacetimes [27j [28] [29] EH] [31] [32] [33] . In view of these reasons, it may be of special interest to investigate 
electromagnetic waves in a plasma in the dS space. 

Recently, Buzzi, et.al. [33J [33], using the 3+1 formulation [35] [35] [37] [35] , described a general rela- 
tivistic version of two-fluid formulation of plasma physics and investigated the nature of plasma waves 
(transverse waves in [33], and longitudinal waves in [34] ) near the horizon of the Schwarzschild black hole. 
In this paper we apply the formalism of Buzzi et.al. |33] to investigate the transverse electromagnetic 
waves propagating in a plasma close to the (cosmological) event horizon of the pure dS space. 

This paper is organized as follows. In section[2]we summarize the 3+1 formulation of general relativity. 
In section[3jwe describe the horizon plasma governing equations. In section[4]we consider one-dimensional 
wave propagation in the radial z (Rindler coordinate system) direction. We linearize the equations in 
section [5] by considering a small perturbation to fields and fluid parameters. In section [6] we discuss 
the local or mean-field approximation for the lapse function a and obtain a dispersion relation for the 
transverse wave. In section [7] we describe a procedure for solving the dispersion relation numerically. In 
sections[8jwe present our results. Finally, in section[9]we give our remarks. We use units G = c = ks = !• 
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2 3+1 Spacetime Formalism 

In de sitter space, the simplest solution for the Einstein held equations with = is written as 

ds 2 = g^dx^dx" 

= -( 1 - r i) dt ' 2 +( 1 - r i) ^ + (1) 

Here, I is the curvature radius of the dS space [i.e., A = -p- is the positive cosmological constant], 
dfl 2 , represents a unit 2-sphere, and the nonangular coordinates range according to < r < £ and 
— oo < t < oo. The boundary at r = t describes a cosmological horizon for an observer located at r = 0. 

An absolute three-dimensional space defined by the hypersurfaces of constant universal time t is 
described by the metric 

ds 2 = g l3 dx l dx j = (l - ^) dr 2 + r 2 d?L 2 2 . (2) 

The indices i, j range over 1, 2, 3 and refer to coordinates in absolute space. The Fiducial Observers 
(FIDOs), the observers remaining at rest with respect to this absolute space, measure their proper time 
r using clocks that they carry with them and make local measurements of physical quantities. Then all 
their measured quantities are defined as FIDO locally measured quantities and all rates measured by 
them are measured using FIDO proper time. The FIDOs use a local Cartesian coordinate system with 
unit basis vectors tangent to the coordinate line 

For a spacetime viewpoint rather than a 3 + 1 split of spacetime, the set of orthonormal vectors also 
includes the basis vector for the time coordinate given by 



did 
6fi ~ lb- ~ adi' 

where a is the lapse function (or redshift factor) defined by 

dr 



(4) 



The gravitational acceleration felt by a FIDO is given by [35, 36, 37, 38 



1 r 

a = -Vina = — -r e^, (6) 
a r 

while the rate of change of any scalar physical quantity or any three-dimensional vector or tensor, as 
measured by a FIDO, is defined by the convective derivative 

D fid 



v being the velocity of a fluid as measured locally by a FIDO. 
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3 Two-Fluid Plasma Equations in 3+1 Formalism 

We consider a two-component plasma consisting of electrons and either positrons or ions. In the 3+1 
notation, the continuity equation for each of the fluid species is 

d 

— (j s n s ) + V • {aj s n s v s ) = 0. (8) 

where s is 1 for electrons and 2 for positrons (or ions). For a perfect relativistic fluid of species s in 

three-dimensions, the energy density e s , the momentum density S s , and stress-energy tensor W% are 
given by 

e s = 7 2 (e s + P s v 2 ), S s - 7s 2 ( £s + P s )v s , Wl k = 7s 2 (e s + P.)«Z«J + P s g^ k . (9) 

where v s is the fluid velocity, n s is the number density, P s is the pressure, and e s is the total energy 
density defined by 

e s = m s n s + P s /{j g - 1). (10) 

The gas constant 7s is 4/3 for T — > oo and 5/3 for T — * 0. 

Using the conservation of entropy, the equation of state can be expressed by 



D ( P. 



0, (11) 



where D/Dt = (l/a)d/dt + v s ■ V. The full equation of state for a relativistic fluid, as measured in the 
fluid's rest frame, is as follows [551141?] : 

P s iH { 2 iy {im s n s /P s )' 



e — m s n s + m s n s 



m s n s iHP(im s n s /P s ) 



(12) 



where the (x) are Hankel functions. 

The quantities of ([9]) in the electromagnetic field are expressed by 

e s = — (E 2 +B 2 ), S s = — E x B, 

8tt ; ' s 4tt 

Wi k = —(F, 2 +B 2 )g jk ~—(E 1 E k + B3B k ). (13) 

87T 47T 

The equations for the conservation of energy and momentum are respectively given by [351 1361 137] 

-JU = -V-S s + 2a-S s , (14) 
a at 

Id 1 ~ 

- — S s = e s a - -V • (qW, ■ (15) 

a at a 

When the two-fluid plasma couples to the electromagnetic fields, the Maxwell's equations take the fol- 
lowing 3+1 form: 

V-B = 0, (16) 

V • E = 4trj, (17) 

— = -Vx(aE), (18) 
d~E 

— = V x (aB) - 47raJ, (19) 



where the charge and current densities are respectively defined by 

a = ^2~fsq s n s , J = ^7 s g s n s v s . (20) 

s s 

Using (|10p and (fT6lfT9|) . the energy and momentum conservation equations (fT4|) and (fT5|) can be 
rewritten for each species s in the form 

a at act 

+ ls q s n s V ■ v s + 2^ 2 {e s + P s )a ■ v s = 0, (21) 

7 2 (e s + P.) r~ + v s • v) v s + VP S - 7,g.n a (E + v s x B) 
+v s ( 7s9s n s E ■ v s + -|-P S J + 7 s 2 (e s + P)[v s (v s ■ a) - a] = 0. (22) 



a dt 

Although these equations are valid in a FIDO frame, they reduce for a — 1 to the corresponding special 
relativistic equations [41] which are valid in a frame in which both fluids are at rest. The transformation 
from the FIDO frame to the comoving (fluid) frame involves a boost velocity, which is simply the freefall 
velocity, given by 

w ff = (l-a 2 )i (23) 

Then the relativistic Lorentz factor 7boost = (1 — W I) -1 ^ 2 = 1/ Q - 

For a good approximation near the horizon, we write the dS metric in the Rindler coordinate system 
as follows: 

ds 2 = - (l-^dt 2 +dx 2 +dy 2 +dz 2 , (24) 

where 

„2\ 1/2 



The standard lapse function in Rindler coordinates becomes a — z/2rh, where r/, = £ is the location of 
the cosmological event horizon. 

4 Radial Wave Propagation in One-Dimension 

We consider one-dimensional wave propagation in the radial z direction and introduce the complex 
variables 

v sz (z,t) = u s (z,t), v s (z,t) = v sx (z,t) +iv sy (z,t), 
B(z,t)=B x (z,t) + iB y (z,t), E(z,t) = E x (z,t) + iE y (z,t). (26) 



Then 



v sx B y -v sy B x = ^(v s B*-v*B), 

v sx E y -v sy E x = ^{v s E* -v* s E), (27) 
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where the * denotes the complex conjugate. The continuity equation |8]) takes the form 

d d 

■7^(l s n s ) + —(aj s n s ii s ) = 0, (28) 



while Poisson's equation (|17[) becomes 



dz 

The and e.g components of (fT8"|) and (fTi?)) give 



47r((7ini7i + 92^272)- (29) 



OE ( 

i— — = —a — a I B - i47reo;(72n2W2 — 7i?iii | i)- (31) 



9i \dz 

Differentiating equation (f3Tj) with respect to £ and using ([30]) . we obtain 

<9 2 3a d d 2 1 \ d 
a2 ^ + ^- h d- Z -W + 4rl) E = ^™d-t {n2l2V2 ~ nillVl) - (32) 

The transverse component of the momentum conservation equation is obtained from the and eg 
components of (|22f as follows: 

Dv s ( 1 dP s \ 

Ps-p^r = q s n s j s (E - iv s B z + \u s B) - u s v s p s a - v s lq s n s j s E, ■ v s + J , (33) 

where 

Ev s = \{Ev* s + E*v s ) + E z u s 



and p s is the total energy density defined by 

Ps = j 2 s (e s + P s ) = 7s 2 (m s n s + T g P s ) (34) 

with T g = 7s/(7 s - !)• 

5 Linearization 

We use perturbation method to linearize the equations derived in the preceding section by introducing 
the quantities 

u s (z,t) = u os (z) + 5u s (z,t), v s (z,t) = 6v s (z,t), 

n s {z,t) = n os (z) + 5n s (z,t), P s (z,t) = P os (z) + 6P s (z,t), 

p s {z,t) = p os (z) + Sp s (z,t), E(z,t) = 5E(z,t), 

B z (z,t) = B (z) + 6B z (z,t), B(z,t)=SB(z,t), (35) 

where magnetic field is chosen to lie along the radial direction. The relativistic Lorentz factor is also 
linearized such that 

7s = Jos + Sis, where j os = (l - u„ s ) 3 , 6j s = ^ s u os ■ 5u s . (36) 
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Near the horizon the unperturbed radial velocity for each species as measured by a FIDO along is 
assumed to be the freefall velocity so that 

u os (z) = v s (z) = [l-a 2 (z)} 1 2. (37) 

It follows, from the continuity equation (|28|) . that 

r 2 a"f os n os u os = const. = r\a h ^ h n h u h , 

where the values with a subscript h are the limiting values at the horizon. The freefall velocity at the 
horizon becomes unity so that Uh = 1- Since u os — vg, 7 os = 1/a; hence, aj os = ah^ih = 1- Also, because 
— r/rh, the number density for each species can be written as follows: 

n os (z) = rihsVfi 3 . (38) 

The equation of state (fTTj) and (|3"8)) lead to write the unperturbed pressure, in terms of the freefall 
velocity, as follows: 

P os (z) = P hs v- 37g . (39) 
Since P os = kBn os T osl then with kg = 1, the temperature profile is 

T os = T hs v s 3 ^- 1] (z). (40) 

The unperturbed magnetic field is purely in the radial direction and it does not experience effects of 
spatial curvature. From the flux conservation V ■ B G = it follows that 

r 2 B (r) = const. 

One can obtain from this the unperturbed magnetic field, in terms of the freefall velocity, as follows: 

B [z) = B h v a 2 , (41) 



where vq = [1 — a 2 (z)] 1 / 2 . Since 



dvff a 1 
dz 2r h vg ' 



(42) 



we have 

duni a 1 dB n a B n 



dz 2rh vs ' dz th Vg ' 

dn os _ dP os _ 3a 7 g F os 

dz 2rh ' dz 2rh vi 

When the linearized variables from (|35[) and (|36p are substituted into the continuity equation and 
products of perturbation terms are neglected, the result gives 

d d u os 2 du os \ r ( d 1 



7 " ( dl + UosU dz + 2^ + ^ a ^z~ ) 5n ° + { a dz + 2^J ^sU os ) 



u — \-g— | 1 g( 1 dn ° s | 3 7 2 u du ° s 
dt dz 2rh \n os dz os dz 



Su s = 0. (44) 



In the similar way, we obtain from the conservation of entropy, (jlip , 



SP S = 1 -^5n s , (45) 

f^OS 



and from the total energy density, (fMj) , 

6 Ps = ^-(l + loslgPoA §ns + 2 Uosl 2 os p os 5u Sl (46) 
where p os = j 2 s (m s n os +T g P os ). Linearizing the transverse part of the momentum conservation equation, 



differentiating it with respect to t and then substituting from (|30|) . we obtain 

d d u os iaq s j 0S n 0S B \ d5v s 



au n 



dz dt 2r/j p os ) dt 



aq s ^ os n os ( d d u os , 



Pos \ dz dt 2r h 
When linearized, Poisson's equation (f29|) and ([32]) respectively give 

= 47re(n o2 7o2 - "-oi7oi) + 47re(7 2<5n 2 - 7 i5ni) 

+47re(n o2 M o2 7o 2 (5u 2 - n ol u ol 'yl l 5ux), (48) 
/ 2 a 2 3a d d 2 1 \ r ( dSv 2 d8v x \ , . 

{ a d^ + ^- h dz-w + ^i) 6E = 4nea { n ° 21o2 ^r ~ nol7ol ^rJ • (49) 

6 Dispersion Relation 

Our consideration effects on a local scale for which the distance from the horizon does not vary signifi- 
cantly. We use a local (or mean-field) approximation for the lapse function and hence for the equilibrium 
fields and fluid quantities. If the plasma is situated relatively close to the horizon, a 2 <g; 1, then a 
relatively small change in distance z will make a significant difference to the magnitude of a. Thus it is 
important to choose a sufficiently small range in z so that a does not vary much. 

We consider thin layers in the direction, each layer with its own a , where a is some mean value of 
a within a particular layer. Then a more complete picture can be built up by considering a large number 
of layers within a chosen range of a values. 

The local approximation imposes the restriction that the wavelength must be smaller in magnitude 
than the scale of the gradient of the lapse function a, i.e., A < (da/dz)^ 1 — 2£, or equivalently, k > (tt/£). 

One of the disadvantages of the hydrodynamical approach is that it is essentially a bulk, fluid approach 
and therefore the microscopic behavior of the two-fluid plasma is treated in a somewhat approximate 
manner via the equation of state. It means that the results are really only strictly valid in the long 
wavelength limit. However, the restriction, imposed by the local approximation, on the wavelength is not 
too severe and permits the consideration of intermediate to long wavelengths so that the small k limit is 
still valid. 

In the local approximation for a, a ~ a a is valid within a particular layer. Hence, the unperturbed 
fields and fluid quantities and their derivatives, which are functions of a, take on their corresponding 
"mean- field" values for a given a Q . Then the coefficients in (|44j) . (|47f and (|48|) are constants within each 
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layer with respect to a (and therefore z as well). Hence, it is possible to Fourier transform the equations 
with respect to z, assuming plane-wave-type solutions for the perturbations of the form ~ e (kz—ut) f or 
each a layer. 

When Fourier transformed, (l47l) and (l49l) turn out to be 



' a o k(a o k~i3/2r h )~L0 2 -l/(2r h ) 2 ' 

( , . iu os a q s "f 0S n 0S B \ . q s 

u a ku os — u+ - 1 ov s - ia a Q ku os - to — - — \SE = 0. (51) 

\ 2r h p os J pos \ 2r h J 

The dispersion relation for the transverse electromagnetic wave modes may be written as 



a 2 J ul 1 {u-u ol K±) ^ uj 2 p2 (lj - u o2 K±) \ 
° \ (u iK T - ui - a ui c i) (u 2K T - u; + a Lu c2 ) \ 



2r h J ~ (2r h ) 

for either the electron-positron or electron-ion plasma. Here, K± = a k ± i/2r^, lu cs — e"f sn sB / p, 



and ujp S = \J '4Tre 2 ^fos n os/ Pos- The cyclotron frequency co cs , as well as the plasma frequency uj ps , is frame 
independent. Although the fluid quantities are measured in the fluid frame, the field B Q is measured in 
the FIDO frame. Hence, the factors of j s do not cancel out explicitly. The transformation B Q — > j sB 
boosts the fluid frame for either fluid and thereby cancels the 7 os factors. The + and — denote the left L 
and right R modes, respectively. The complex conjugate of the dispersion relation for the R mode gives 
the dispersion relation for the L mode. In the special relativistic case, the two modes have the same 
dispersion relation. 

7 Numerical Solution Modes 



The dispersion relations (|52|) are complicated enough even in the simplest cases for the electron-positron 
plasma where both species are assumed to have the same equilibrium parameters, and an analytical 
solution is cumbersome and unprofitable. We therefore solve numerically the dispersion relation in order 
to determine all the physically meaningful modes for the transverse waves. We put the equations in the 
form of a matrix equation as follows: 

(A-kI)X = 0, (53) 

where the eigenvalue is chosen to be the wave number k, the eigenvector X is given by the relevant set 
of perturbations, and / is the identity matrix. 

In order to write the perturbation equations in an appropriate form, we introduce the following set 
of dimensionless variables: 

kh 





U) 




kc 


Ul = 




k 


? 




_ Su s 




_ 5v s 


Su s - 


~ ? 


v s 






SB 




SE 


SB 




E 


~B~o 



2r h u>* 
Sn. 



Sh s = 

SE Z = 6 -^. (54) 

For an electron-positron plasma, oj p \ — uj P 2 and uj c \ = co C 2', so, a->* is defined as 

Alfven modes, 

co* = ^ (55) 
(2lu 2 + lj 2 )i , high frequency modes, 
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where u> p — ^/uj p iujp2 and u c — \Ju) c \<lo c i. However, for the case of an electron-ion plasma, the plasma 
frequency and the cyclotron frequency are different for each fluid; so, the choice of is a more complicated 
matter. We assume, for simplicity, that 



^(k>ci + u>&) * , Alfven modes, 

(w^ + uj^ 2 )K high frequency modes, 



where ux* a = (2w| a + ui 2 cs ). 

The dimensionless eigenvector for the transverse set of equations is 



X 



Lraiisvcrsc 



Svi 

SB 
SE 



When linearized and Fourier transformed, equations (|30p and (|3ip turn out to be 

k 



1 



SE+—6B = 0, 
2r h a j a a 



™6E= (k--L- 
a \ 2rha Q 



SB + 47re(7 o2 n o2 (5w2 - loin i5vi). 



Using ([53]). we write (JSJ), 



and (|59|) in the dimensionless form: 



fccH;., = I — 

, u ns V e / u„ s oj 



a a I \ e / 



— <W , 



(!) 



fe<5£ = 
fc<5£> = u i 



a„ 



■OVi — U 2 0V2 H SB + VJJOhi. 



U) c \LO* U) C 2U}„ 

These are the equations in the required form to be used as input to ([5 



(56) 



(57) 



(58) 
(59) 

(60) 
(61) 
(62) 



8 Results 

We carried out the numerical analysis using the well known MATLAB. We have considered both the 
electron-positron plasma and the electron-ion plasma. The limiting horizon values for the electron- 
positron plasma are taken to be 

4 



n hs = 10 18 cm^, T hs = 10 iu K, B h = 3 x 10 b G, and lg 



,10 



(63) 



For the electron-ion plasma, the ions are essentially non relativistic, and the limiting horizon values are 
chosen to be 

n hl = 10 18 cm- 3 , T hl = 10 10 K; n h2 = 10 15 cm~ 3 , T h2 = 10 12 K. (64) 

The equilibrium magnetic field has the same value as it has for the electron-positron case. The gas 
constant is "f g = 4/3. 
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8.1 Alfven Modes 

8.1.1 Electron-Positron Plasma 

For the ultrarelativistic electron-positron plasma in the special relativistic case, only one purely real 
Alfven mode exists |41j , while for the Schwarzschild case there are two Alfven modes [33] . In our analysis 
we find three Alfven modes, as shown in Fig. 1 and Fig. 2., for the electron-positron plasma. The Alfven 
modes in de Sitter space are interesting in that, there exists three Alfven modes for the electron-positron 
plasma compared with four modes for the electron-ion plasma. These three modes for the electron- 
positron plasma coalesce into a single mode on taking the special relativistic limit, giving the result of 
ref. [41]. Since we have used the convention e lkz = e 1 l Rc ( k )+ lIm ( k )} ; the damping corresponds to Im(fc) > 
and growth to Im(fc) < 0. 

8.1.2 Electron-Ion Plasma 

In this case four modes are found, two of which are growth and the other two are damped. The modes 
shown in Fig. 3 and Fig. 4 are damped and the remaining two modes shown in Fig. 5 and Fig. 6 are 
growth. The first two modes arc equivalent to the modes shown in Fig. 1 and the other two are equivalent 
to the modes shown in Fig. 2 for electron-positron plasma. 

The differences in the magnitudes of the cj c i and oj C 2 for the first two modes apparently lead to take 
the frequencies from their negative (and therefore unphysical) values for the electron-positron case to 
positive physical values for the electron-ion case. These changes are thus because of the difference in 
mass and density factors as between the positrons and ions. 

These four modes for electron-ion plasma are equivalent to those of the Schwarzschild case [33]. It is 
evident that the growth and damping rates are independent of the frequency, but depended only on the 
value of a . 

8.2 High Frequency Modes 
8.2.1 Electron-Positron Plasma 

In this case three high frequency electromagnetic modes are found for the electron-positron plasma, as 
shown in Figs. 7-9. High frequency modes in the horizon of dS space for each fluid are interesting in 
that all the modes are both damping and growth modes. 

The two modes, shown in Figs. 7 and 8, are similar and both modes are damping and growth very 
near to the horizon. These two modes are growth for most of the frequency domain but shows damped 
for lower frequencies as uj — ► .04 and a — ► 0. Thus at a distance from the horizon corresponding to 
a — > 0.2 it appears that energy is no longer fed into wave mode by the gravitational field but begins to 
be drained from the waves. The third mode shown in Fig. 9 is also growth and damping mode. This 
mode is damped for most of the frequency domain but shows growth for lower frequencies as uj — ► .04 
and a — ► 0. These three modes are equivalent to the three modes of ref. [33] for this case. Also these 
three modes coalesce with a single modes in the special relativistic case [H] as a — > 1. 
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8.2.2 Electron-Ion Plasma 

Similar as for the electron-positron plasma, the electron-ion plasma has three high frequency modes. Two 
of these, shown in Fig. 10 and Fig. 11, shows damping and growth very near to the horizon for lower 
frequencies, and for upper frequencies they are growth modes. These two modes are growth for higher 
frequency but shows damped for lower frequencies as to — > .02 and a Q — > 0. The third mode, shown in 
Fig. 12, is growth for lower frequencies as u> — ► .02 and a — ► and shows damped for higher frequencies. 
These three modes are similar to the modes for high frequency electron-ion plasma in the Schwarzschild 
case [33] . 

9 Concluding Remarks 

The main concern of this study has been exclusively the investigation, within the local approximation, 
of Alfven and high frequency transverse electromagnetic waves in a two-plasma in the purely de Sitter 
space. We derive the dispersion relations for the Alfven and high frequency electromagnetic waves by 
using a local approximation and give their numerical solutions. In the limit £ — > our results reduce 
to that in special relativity as obtained by Sakai and Kawata [41] (i.e., only one purely real mode for 
Alfven and high frequency electromagnetic waves). In contrast to the work of Sakai and Kawata [41) . new 
modes (damped or growth) arise for the Alfven and high frequency electromagnetic waves in the pure dS 
space. In our work all the modes for Alfven waves are either damped or growing, but for high frequency 
electromagnetic waves all the modes are both damped and growing. This is because of the singularity of 
the de Sitter space. For the electron-positron plasma, the damping and growth rates are similar with the 
electron-ion plasma but different by several orders of magnitude, compared with the real components of 
the wave number. For both the fluid components the damping and growth rates are obviously frequency 
independent, but are dependent on the radial distance from the horizon as denoted by the mean value 
of the lapse function a . This is of course not for the case of the high frequency waves. In that case 
the rate of damping or growth is dependent on both frequency and radial distance from the horizon. 
Damped modes demonstrate, at least in this approximation, that energy is being drained from some of 
the waves by the gravitational field. The majority of the modes are growth rates and that indicate that 
the gravitational field is feeding energy into the waves. 

In the light of recent astronomical observations, it has been suggested that our universe will asymp- 
totically approach a de Sitter space pQ. Hence, aspects of the de Sitter space might be of interest in a 
broader context. Our study of plasmas in the de Sitter space is thus well motivated. 
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Figure 1: Top: Real part of Alfven mode for the electron-positron plasma. Middle: Imaginary part of 
Alfven damped mode. Bottom: Imaginary part of Alfven growth mode. 
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Figure 2: Left: Real part of Alfven mode for the electron-positron plasma. Right: Imaginary part of 
Alfven growth mode. 




Figure 3: Left: Real part of Alfven mode for the electron-ion plasma. Right: Imaginary part of Alfven 
damped mode. 




Figure 4: Left: Real part of Alfven mode for the electron-ion plasma. Right: Imaginary part of Alfven 
damped mode. 
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Figure 5: Left: Real part of Alfven mode for the electron-ion plasma. Right: Imaginary part of Alfven 
growth mode. 




Figure 6: Left: Real part of Alfven mode for the electron-ion plasma. Right: Imaginary part of Alfven 
growth mode. 




Figure 7: Left: Real part of high frequency mode for the electron-positron plasma. Right: Imaginary part 
of high frequency damping and growth mode. 
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Figure 8: Left: Real part of high frequency mode for the electron-positron plasma. Right: Imaginary part 
of high frequency damping and growth mode. 




Figure 9: Left: Real part of high frequency mode for the electron-positron plasma. Right: Imaginary part 
of high frequency growth and damping mode. 




Figure 10: Left: Real part of high frequency mode for the electron-ion plasma. Right: Imaginary part of 
high frequency damping and growth mode. 
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Figure 11: Left: Real part of high frequency mode for the electron-ion plasma. Right: Imaginary part of 
high frequency damping and growth mode. 




Figure 12: Left: Real part of high frequency mode for the electron-ion plasma. Right: Imaginary part of 
high frequency growth and damping mode. 
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